Hundred photon microwave ionization of Rydberg atoms in a static electric field 
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We present analytical and numerical results for the microwave excitation of nonhydrogenic atoms 
in a static electric field when up to 1000 photons are required to ionize an atom. For small microwave 
fields, dynamical localization in photon number leads to exponentially small ionization while above 
quantum derealization border ionization goes in a diffusive way. For alkali atoms in a static field 
the ionization border is much lower than in hydrogen due to internal chaos. 
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After two decades of investigations initiated by the pi- 
oneer experiment of Bayfield and Koch in 1974 jl), the 
main features of the microwave ionization of highly ex- 
cited hydrogen atoms are well understood Gener- 
ally, this ionization takes place due to emergence of chaos, 
in the corresponding classical problem, above a certain 
microwave intensity threshold. In the chaotic regime, 
quantum excitation proceeds via a diffusive absorption 
and reemission of field photons, which eventually leads 
to ionization. For a real atom the quantum process can 
be either close to the classical diffusion, or strongly sup- 
pressed when the field e is below a quantum derealiza- 
tion border i q . This suppression is due to dynamical 
localization of classical chaos which is produced by quan- 
tum interference. In the localized case the distribution 
over the levels drops exponentially with the photon num- 
ber. This dynamical localization is analogous to the An- 
derson localization in quasi one-dimensional disordered 
solids, with the photon number playing the role of lat- 
tice index. The quantum derealization takes place if 
the localization length > Nj, where, in atomic units, 
Nj = l/(2ngw) is the number of photons required for 
ionization, with no being the initial principal quantum 
number and lj the field frequency. 

An absolutely different scenario for the microwave ion- 
ization of alkali Rydberg atoms had been proposed by 
Gallagher et al. based on experimental results in the 
regime with small rescaled frequency loq — louq -C 1 |||| . 
According to this scenario, ionization appears due to a 
chain of consecutive Landau-Zener transitions between 
nearby levels, which eventually brings the electron into 
the continuum. In the presence of quantum defects the 
Stark manifolds exhibit a structure of avoided crossings 
only for sufficiently strong fields with e > 1/t1q. Indeed, 
the experiments (5]|(| were in agreement with the e ~ 
1 /tiq dependence for the microwave ionization threshold 
and not with the static field border e s = l/(9ng)- This 
was considered as an experimental confirmation of the 
above scenario. However, this picture doesn't explain in 



fact how the propagation actually occurs via the chain of 
these transitions and why the overlapping of two nearby 
levels guarantees that the electron will pass through the 
whole chain. The comparison between the thresholds for 
hydrogen jj] and for alkali atoms 1^-0 clearly shows that 
the border is lower in the latter case and therefore it is 
related to the quantum defects Si of alkali atoms. Since 
these defects are different from zero only for orbital mo- 
mentum I < 3, the situation is purely quantum and can- 
not be treated by the quasi-classical approach used for 
hydrogen atom M . This is the reason why so long after 
experiments have been made |5|,^j no detailed theory has 
been developed. To demonstrate a theoretical difficulty 
we note that hundreds of photons are required to ionize 
atoms in this regime (Ni = no/2ujQ = 300 for no = 60 
and ujq = 0.1). 

In this Letter we propose another mechanism for the 
microwave ionization of Rydberg atoms in a static electric 
field which is qualitatively different from the Gallagher 
et al. scenario. We argue that ionization in this case 
is not due to Landau-Zener transitions but to a quan- 
tum diffusive excitation in energy E (or in the photon 
number N = E/lu). Such quantum diffusion in the low 
frequency regime 1/uq < ujq <C 1 becomes possible due 
to appearence of quantum chaos for Rydberg atoms in 
a static electric field. Indeed, a recent theoretical study 
H showed that the level spacing statistics P(s) in such 
atoms, for a sufficiently strong static field, is described by 
the Random Matrix Theory. These results indicate also 
a chaotic structure of cigenstates. This internal quan- 
tum chaos can lead to a diffusive excitation in energy 
even for a quite weak microwave field. In this respect, 
the situation is different from the hydrogen atom where 
chaos could appear only above some classical field thresh- 
old. Another consequence of internal chaos is an increase 
of density of effectively coupled states p c . This gives 
a larger localization length £ q (£ q oc p c ) and therefore 
significantly decreases the derealization border as com- 
pared to the hydrogenic case. A similar effect of chaotic 
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enhancement of localization length due to internal chaos 
has been studied recently for hydrogen atoms in magnetic 
and microwave fields |@. 




with the core, which are responsible for diffusion, hap- 
pen rarely, with a frequency of classical precession in I, 
which is lo s q — 3e s o <C 1. Instead, for ujq > 1 the ratio 
D q /D w 1/3cl>q^ 3 is similar to the hydrogen case |pl| . 
Notice that the dependence of D q /D on ujo is very fiat in 
the interval 0.02 < ujq < 0.3. The origin of this fact is not 
quite clear since asymptotically ,for small u)q, one should 
expect D q oc uj 2 |§,|l(j. (Apparently, this plateau in fre- 
quency appears due to one-photon transitions between 
levels in the Stark multiplet. As a result its rescaled size 
is expected to be lo s q w 3e s o, which is in approximate 
agreement with the data of Fig.l.) We also checked that 
a change of uq from 60 to 28 at the same classically scaled 
parameters did not affect the ratio D q /D$. 



FIG. 1. Diffusion rate ratio D q /Do as a function of mi- 
crowave frequency u>o for Rb (circles), Na (triangles) and Li 
(squares) at no = 60, e s o = 0.02. The full line is drawn to 
guide the eye. The upper insert gives < (AN4,) 2 > vs. r at 
eo = 0.03, u>o = 0.1: diffusive excitation for Rb (full curve 
D q /Do = 0.005) and small oscillations for hydrogen (dot- 
ted curve with 100 times magnification). The lower insert 
illustrates the weak dependence of D q /Do on e s o for Rb at 
cjo ~ 0.1, eo > 0.015 (crosses). 

To check the above picture of quantum photonic dif- 
fusion, we numerically studied the excitation of alkali 
Rydberg atoms (Rb, Na, Li) in a static electric field e s 
and parallel, linearly polarized, microwave field esinwt, 
for magnetic quantum number m = 0. We chose the 
rescaled value e s o = e s n,Q as 0.02 so that the statis- 
tics P(s) for levels with 55 < no < 72 in Rb and 
Na was close to the RMT results. At the same time 
P(s) for Li was closer to Poisson statistics due to a 
smaller value of quantum defects fll|| . The investiga- 
tion of time evolution, in the eigenbasis of the unper- 
turbed problem (e = 0), showed that an initial eigen- 
state with energy E spreads diffusively over the unper- 
turbed energies E\, namely the square variance of the 
photon number a =< (AA^,) 2 >=< (Ex - E ) 2 > /uj 2 
initially grows linearly with the number of microwave pe- 
riods r. Our quantum simulation allows to determine the 
value of quantum diffusion rate in energy per unit time 
D q =< (AE) 2 > /At. This rate D q can be compared 
with the diffusion rate in hydrogen at uiq = 1, given by 
Dq = e 2 no/2 Q. The dependence of the ratio D q /Do on 
the frequency w , for rescaled fields 0.005 < e < 0.03 
(eo = en o)i is shown in Fig. 1. 

In our computations we consider the most interesting 
frequency region ujq < 0.5 where ionization in hydro- 
gen is close to the static field ionization. In our opinion 
the ratio D q /Do in Fig. 1 is small because collisions 




-60 -40 -20 O 20 40 60 SO 

N<4 



FIG. 2. Probability distribution over the eigenstates at 
e = (full line) and in one-photon intervals (circles) 
for Rb: n = 60, e s0 = 0.02, u = 0.1, eo = 0.005, 
180 < r < 200, A,/ A) = 0.0051,4, = 9.1. The 
straight line shows the fit for exponential localization with 
£ q N ~ 13.3. The one-photon probabilities are also shown, in 
the same conditions, for Na (triangles, shifted down by 10 3 , 
£ q = 4.8, £ q N — 6.4) and Li (squares, shifted down by 10 6 , 
4 = 3.1,4jv = 5.4). 

The above diffusive excitation in energy induced by 
internal chaos may eventually be localized at long times 
due to quantum interference effects in a way similar to 
photonic localization in a complex molecular spectrum 
Jl2|,P|. The localization length £ q , expressed in the num- 
ber of photons, is proportional to the one-photon tran- 
sition rate T and to the density of coupled states p c : 
£ q ~Tp c pi. Since T - D q /uj 2 we obtain 
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Here the length . 
in hydrogen - 



is expressed via the localization length 
3.3eo?T.o at wo = l. Notice that in the 
non-hydrogenic case the density of states is p c = Uq due 
to internal chaos while in hydrogen the effective density 
is smaller, that is p — due to existence of an addi- 
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tional integral of motion |2j]. The above expression for 
l q is valid for £ q > 1 and uip c = uqUq > 1. The local- 
ization leads to an exponential decay of probability dis- 
tribution IV'aH 2 ~ exp( — 2|A^|/£q) in the photon number 
N^ = (E x -E )/lj. 
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FIG. 3. The numerically computed localization length £ q M 
vs. theoretical £ q (Eq. (1)) for Rb (points) and Na (triangles) 
at n = 60, e s0 = 0.02 and 0.02 < u < 0.5, 0.003 < e < 0.01. 
The straight line gives £ q N = £ q , while the numerical average 
is < £ qN /£ q >= 1.17 ±0.28. 

In order to check the theoretical prediction (1) we nu- 
merically computed the quantum evolution following it 
up to 200 microwave periods. The probability distri- 
bution f\ over the eigenstates of the static field prob- 
lem (e = 0) with energies E\, was averaged over 10 - 
20 periods to suppress the fluctuations. An initial state 
at e s o = 0.02 was chosen as an eigenstate with energy 
E\ a k Eq = l/2ng and no = 60. The system param- 
eters were varied in the intervals: 0.02 < uiq < 0.5, 
0.003 < e < 0.03 and 60 < A/ < 1500. The total 
basis included up to 1150 states. 

A typical example of stationary distribution / is shown 
in Fig. 2. It clearly demonstrates exponential localiza- 
tion of diffusive excitation. One can note that among 
the three cases shown (Rb, Na, Li) the most localized is 
the case of Li, for which the quantum defect is minimal 
and therefore the internal chaos is the most weak. Notice 
also that in the case of Na and Li, for Jn < 10~ 7 , the 
probability /at starts to decay in a much slower way with 
l q » 25 (Fig. 2). We attribute this effect to a significant 
modification of hydrogenic basis on highly excited levels 
where a static field becomes quite strong and tunneling 
effects for probability decay should be taken into account. 

In order to find the localization length £ q we first com- 
pute the total probabilities /jy = |^Af| 2 in one-photon 
intervals [N^ — 1/2, + 1/2] around integer values of 
— E\/oj and then extract the numerical £ q N value 
from the least square fit for In /at. The analysis of the 



numerical data £ q N shown in Fig. 3 confirms the theoret- 
ical prediction (1) for £ q with D q rates taken from Fig.l. 

Equation (1) allows to determine the quantum delo- 
calization border above which localization effects become 
unimportant and ionization goes in a diffusive way. This 
happens for £ q > A/, which gives the quantum delocal- 
ization border for the rescaled field e q = e q n^: 

e q = QA^/^y^Do/Dq. (2) 

Our numerical data indeed show that above this border 
complete dclocalization takes place, contrarily to the case 
of hydrogen atom at the same field parameters (Fig. 4). 
The strong fluctuactions in the hydrogenic distribution 
/at indicate that it is quite inhomogeneous inside the 
atomic shells. 
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FIG. 4. Same as Fig. 2 but in a delocalized case for Rb, with 
e = 0.03 > e q = 0.028 and 40 < r < 50. Full triangles give 
one-photon interval probability distribution }m for hydrogen 
atom at the same system parameters (dashed lines show the 
distribution envelope). 

For fixed luq the delocalization border (2) scales as 
iq ~ 1/^q, namely it is ^/hq times smaller than for 
hydrogen at u>o ~ 1. This drop of e q is related to the 
appearence of internal chaos in nonhydrogenic atoms in 
a static electric field which effectively enhances the in- 
teraction with the microwave radiation. According to 
(2) the ratio of e q to the ionization border proposed 
by Gallagher et al. tQ sa 1/n-o depends only on luq: 
£q/tG = OAyJuoDo/Dq. Our data from Fig.l indicate 
that this ratio varies rather weakly with ujq in the inter- 
val 0.02 < o^o < 0.5 (see Fig. 5). In spite of the fact that 
in this range e q ~ e<3, the physical interpretation is rather 
different from Gallagher et al. scenario. Indeed, the sit- 
uation for e < e q ~ V n o is strongly nonperturbative 
since for e g (2w /no) 1 ^ 2 < e o many photons are absorbed 
with l q > 1. 

The theory developed above allows to understand ion- 
ization of atoms in a static electric field. It is possible to 
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expect that for luq <C 1 the situation will remain similar 
even without static field since its role will be played by a 
slowly varying microwave field. However, one should be 
careful in extending the theory to a zero static field case. 
Indeed, our numerical results show that there, the prob- 
ability distribution in orbital momentum is qualitatively 
different. For example, at eq ~ 0.02, ljq ~ 0.2 only few 
Z-states are mixed while with additional e s o ~ 0.02 prob- 
ability spreads over all accessible I. The localization in 
I space at e s q = had been also discussed for Rb atoms 
in |l3| . The physical reason of this difference is related 
to the fact that the condition luq <C 1 is not sufficient 
to treat the microwave field eo as quasi-static. For that 
one should require ujq <C u> s o ~ 3eo since the precession 
frequency cu s0 determines oscillations in I Q. However, 
this condition is not compatible with the requirement 
loq > 1 /no for the considered region of no. For e s o = 
chaos is induced by the microwave field, the distribution 
in I is nonhomogeneous and a detailed theory for this case 
becomes more complicated as compared to the nonzero 
static field where internal chaos is already present. 



liminary numerical data |llj show that there is a change 
of slope in the probability decay similar to one in Fig. 2 
{l qN « 4 for f N > IQ- 5 and l qN « 25 for f N < 10~ 5 ). 
This means that in this case tunneling ionization plays 
a dominant role while the microwave only slightly in- 
creases its rate. An increase of the principal quantum 
number up to no = 60 strongly suppresses the tunneling 
in a static field and the diffusive microwave excitation 
becomes dominant (Figs. 2,4). The effects of tunneling 
for no ~ 30 will become less important for experiments 
with a shorter interaction time (r ~ 100) on which the 
dynamical localization dominates. It is also possible that 
for such long times as in [[| some noise in a microwave sig- 
nal can strongly affect ionization. The above discussion 
shows that more detailed experimental investigations of 
microwave ionization in a static electric field are highly 
desirable. They will allow to make a detailed test of 
dynamical localization theory in the regime when up to 
1000 photons are required to ionize an atom. 
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FIG. 5. The ratio of the quantum derealization border e 9 
to the border eo ~ 1/nn vs. uiq for Rb (circles) obtained from 
Eq. (2) and data of Fig.l. The full line is drawn to guide the 
eye, crosses refer to the cases of Figs. 2,4. 

Another difficulty for direct comparison with the ex- 
periments is that the latter were mainly done in 
the regime luq < 1/no- There is only one case for Na at 
no = 28, wo.no = 0.76 and e s o « 0.024 which is not far 
from our regime (Fig. 2d in 0). Here the experiment 
gives the ionization border eoex ~ 0.002. This value is 
about 20 times smaller than the quantum derealization 
border given by Eq. (2) with D q /D = 0.0027. We re- 
late this difference to a tunneling ionization in a static 
field which becomes important for very long interaction 
times (3 x 10 5 Kepler periods in ||). Indeed, our pre- 
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